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HIGHLIGHTS 


•  Microstructure  of  a  bi-layer  Ni— YSZ  cermet  is  analysed  from  3D  reconstructions. 

•  Functional  layer  and  substrate  are  quantified  with  statistically  representative  volumes. 

•  Microstructure  properties  of  all  phases  are  determined  for  the  Ni-YSZ  bi-layer. 

•  Extent  of  the  electrochemical  reaction  in  the  support  is  discussed. 

•  Electrode  activation  energy  is  computed  as  function  of  microstructure  characteristics. 
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Microstructure  of  a  cermet  Ni-YSZ  bi-layer  is  analysed  on  the  basis  of  three  dimensional  reconstructions 
obtained  on  both  functional  layer  and  cell  support.  Microstructural  parameters  of  gas,  ionic  and  elec¬ 
tronic  phases  are  determined  in  terms  of  phase  connectivity,  mean  particles  diameter,  particles  size 
distribution,  specific  surface  area,  tortuosity  factor  and  density  of  TPBls.  Microstructural  properties  are 
introduced  in  an  SOEC  cathode  micro  model  that  takes  into  account  the  specific  configuration  of  the  Ni 
—YSZ  composite  bi-layer.  Simulations  show  that  the  extent  of  the  electrochemical  reaction  in  the  support 
is  very  limited.  Moreover,  it  is  found  that  electrode  apparent  activation  energy  is  a  combination  of 
effective  ionic  conduction  and  charge  transfer  in  the  active  functional  layer. 
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1.  Introduction 

Because  of  its  great  potential,  hydrogen  production  by  High 
Temperature  Steam  Electrolysis  (HTSE)  has  received  an  increasing 
interest  in  recent  years.  Hydrogen  produced  by  water  electrolysis 
taking  advantage  of  sustainable  energy  sources  will  allow  limiting 
greenhouse  gas  emission.  However,  this  scheme  will  be  achievable 
only  if  the  efficiency,  durability  and  reliability  of  the  electrolyser  are 
improved. 

The  active  parts  of  an  electrolyser  correspond  to  a  stack  of 
elementary  cells,  named  Solid  Oxide  Electrolysis  Cells  (SOECs). 
SOECs  are  constituted  of  a  three-layered  structure  made  of  ceramic 
materials:  one  dense  electrolyte  ofYttria  Stabilized  Zirconia  (YSZ)  is 
sandwiched  between  two  porous  electrodes.  The  cathode  material 
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is  usually  a  composite  of  Nickel  and  YSZ  (Ni-YSZ).  The  Cathode 
Supported  Cell  (CSC)  configuration  is  composed  of  a  thick  highly 
porous  Ni-YSZ  substrate  on  which  a  thin  electrolyte  and  the 
functional  electrodes  are  layered.  This  geometry  has  been  designed 
to  improve  the  cell  electrochemical  efficiency  by  limiting  the  ionic 
losses  in  the  thin  electrolyte.  Nowadays,  the  planar  CSC  geometry 
allows  achieving  the  best  electrochemical  performance  at  usual 
SOEC  operating  temperatures  (-800  °C)  [1-3]  and  is  generally 
employed  in  stack  preindustrial  demonstrators  [2-4]. 

In  CSC  configuration,  the  cathode  substrate  exhibits  a  coarse 
microstructure  which  is  crucial  for  the  gas  transportation.  Indeed, 
the  thick  substrate  can  potentially  hinder  the  gas  diffusion  and  can 
induce  high  concentration  overpotential  at  high  current  density 
[5-7].  Conversely,  the  functional  layer  presents  a  thinner  micro¬ 
structure  which  must  be  optimized  to  facilitate  the  electrochemical 
process  and  minimize  the  activation  overpotential.  For  this  pur¬ 
pose,  the  density  of  Triple-Phase  Boundary  lengths  (TPBls)  should 
be  increased  as  far  as  possible  since  it  corresponds  to  the  active 
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Fig.  1.  Three  dimensional  rendering  volumes  of  studied  Ni-YSZ  cermet:  (a)  reconstruction  of  36  x  36  x  36  |im3  taken  in  the  middle  of  the  cell  support,  (b)  Reconstruction  of 
27.5  x  27.5  x  27.5  pm3  including  the  functional  layer,  a  part  of  the  cell  support  and  a  part  of  the  electrolyte. 


sites  where  the  electrochemical  reactions  take  place  [8].  In  addi¬ 
tion,  microstructure  of  the  functional  layer  must  be  tailored  to 
promote  the  effective  ionic  conduction  of  the  YSZ  structure.  This 
requirement  implies  to  lower  the  YSZ  tortuosity  factor  so  that  the 
oxygen  ions  conductivity  is  enhanced  in  the  functional  layer. 
Indeed,  it  has  been  shown  that  this  parameter  governs  the  extent  of 
the  electrochemical  reaction  from  the  electrolyte  interface  within 
the  functional  layer  [8]. 

Therefore,  microstructure  properties  of  both  functional  layer 
and  substrate  have  to  be  quantified  precisely  and  correlation  with 
electrode  polarisation  must  be  established.  For  this  purpose,  there 
is  a  strong  interest  to  use  micro  models  able  to  analyse  the  impact 
of  electrode  microstructure  properties  on  activation  and  concen¬ 
tration  overpotentials.  In  a  previous  paper  [8],  the  authors  estab¬ 
lished  such  an  electrochemical  micro  model  and  performed 
parametric  studies  to  evaluate  some  degradation  mechanisms. 

In  this  paper,  it  is  proposed  to  compute  the  microstructure 
features  of  a  typical  CSC  Ni-YSZ  electrode  from  three  dimensional 
reconstructions  obtained  by  X-ray  nano-tomography  [9].  Both 
functional  layer  and  support  are  analysed  and  the  microstructure 
properties  of  gas,  ionic  and  electronic  phases  are  determined  for 
the  two  layers.  For  the  first  time,  the  main  microstructural  char¬ 
acteristics  of  a  Ni-YSZ  cermet  bi-layer  are  quantified  with  a  large 
field  of  view  (~50  pm),  in  terms  of  (i)  tortuosity  factors  of  each 
phase,  (ii)  density  of  TPBls,  (iii)  Particles  Size  Distribution  (PSD)  and 
(iv)  mean  particles  diameters.  The  local  in-house  micro-model  [8]  is 
then  used  to  relate  the  microstructure  properties  to  the  electrode 
polarization.  Impact  of  the  microstructure  on  the  global  electrode 


response  is  analysed.  In  this  frame,  a  special  attention  has  been  paid 
to  estimate  the  possible  electrochemical  activity  of  the  Ni-YSZ 
substrate.  For  this  purpose,  the  extent  of  the  electrochemical  re¬ 
action  from  the  functional  layer  into  the  support  is  calculated  as  a 
function  of  its  real  microstructure  properties.  The  extraction  of  the 
electrode  exchange  density  as  a  function  of  temperature  is  also 
addressed  in  this  article.  It  is  shown  that  the  combination  of  ionic 
conduction  and  charge  transfer  all  along  the  volume  results  in  an 
apparent  activation  energy  or  the  electrode  activation  consistent 
with  the  published  experimental  data. 

2.  Material  and  methods 

2.1.  Sample  and  X-ray  nano -tomography  characterisation 

Sample  investigated  in  the  present  work  was  taken  from  com¬ 
mercial  CSC  supplied  by  the  Forschungszentrum  Jiilich  (FZJ) 
research  centre  [10].  Cathode  is  constituted  by  a  bi-layer  made  of 
50  wt%  Ni  and  50  wt%  8YSZ  (Zr02  stabilised  with  8  mol%  Y2O3).  It 
includes  a  thin  Cathode  Functional  Layer  (CFL)  deposited  onto  a 
thick  Cathode  Current  collector  (CCC).  The  CCC  substrate  is  ~  1  mm 
thick  and  is  used  as  the  cell  structural  support.  The  CFL  and  the 
8YSZ  electrolyte  are  co-sintered  with  the  cell  support  and  present  a 
thickness  of  ~10  pm. 

Synchrotron  X-ray  nano-tomography  experiments  have  been 
carried  out  at  the  European  Synchrotron  Radiation  Facility  (ESRF 
beam  line  ID22)  to  characterize  the  three  dimensional  structure  of 
CFL  and  CCC  layers.  The  tomography  procedure  is  based  on  an 
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optical  phase  contrast  imaging  technique  that  follows  a  hol- 
otomographic  scheme  [11  .  The  method  is  well  adapted  to  separate 
the  porosity  [12]  as  well  as  the  different  solid  phases  having  the 
same  X-ray  absorption  such  as  Ni  and  YSZ  [9].  The  3D  re¬ 
constructions  present  a  large  field  of  view  ( ~  50  pm)  with  a  voxel 
size  of  25  nm.  Further  details  on  data  acquisition,  algorithm  for  3D 
reconstructions  and  image  segmentation  are  reported  in  Villanova 
et  al.  [9].  The  reconstructed  volume  of  the  studied  CFL  and  CCC 
layers  is  also  provided  in  Ref.  [9].  Fig.  la  and  b  illustrates  3D 
rendering  volumes  related  to  the  Ni-YSZ  support  and  functional 
layers.  It  can  be  noticed  that  these  volumes  are  sufficiently  large  to 
be  statistically  representative  for  the  studied  heterogeneous  media. 
Indeed,  the  Representative  Volume  Elements  (RVEs)  are  roughly 
equal  to  35  x  35  x  35  pm3  [12]  and  13  x  13  x  13  pm3  [9  for  the 
current  collector  and  functional  layer  respectively. 

2.2.  Procedures  for  microstructures  quantification 


Complementary  phase 


A  set  of  programs  has  been  implemented  in  Matlab®  software  to 
compute  all  the  microstructural  properties  required  for  the  elec¬ 
trode  modelling.  Methods  and  algorithms  to  determine  the  volume 
fraction,  phase  connectivity  and  “apparent”  tortuosity  factors  have 
been  already  presented  in  Laurencin  et  al.  [12  .A  short  abstract  is 
given  hereafter  while  methods  for  determining  surface  specific 
area,  mean  particle  diameter,  particle  size  distribution  and  length  of 
TPBs  are  described  more  precisely. 

2.2.1.  Phase  connectivity 

For  determining  the  phase  contiguity,  it  is  reminded  that  two 
cubic  voxels  of  the  reconstruction  are  assumed  to  be  connected 
only  if  they  share  a  common  face.  Inversely,  two  neighbouring  cubic 
voxels  are  supposed  to  be  disconnected  if  they  are  only  linked  by  an 
edge  or  a  vertex. 


Fig.  2.  Schematic  representation  of  the  analysed  phase  for  PSD  calculation.  The  volume 
is  covered  by  a  set  of  overlapping  spheres  (which  must  be  tangential  to  the  comple¬ 
mentary  phase).  The  voxels  are  attributed  to  the  largest  sphere. 


2.2.3.  Mean  particle  diameter  and  surface  specific  area 

The  mean  particle  diameter  and  surface  specific  area  can  be 
computed  from  basic  morphological  tools  [17,18].  The  geometrical 
covariogram  I<(h)  can  be  presented  as  a  measure  of  the  intersection 
between  the  phase  and  its  image  after  a  translation  h  of  the  me¬ 
dium.  Considering  the  x-direction,  K(h)  is  expressed  as  followed: 

K(h)  =  fffk(x,y,z)  xk((x+h),y,z)xdV  and 

V(h)  (2) 
lc(xvz)  —  |1  if  Xiongs  to  the  studied  phase  v  ; 

\  -y-  '  if  (x,y,z)  isnotincludedin the studiedphase 


2.2.2.  Tortuosity  factor 

The  “apparent”  tortuosity  factors  of  gas,  ionic  and  electronic 
conducting  phases  are  determined  by  using  a  homogenization 
technique.  For  this  purpose,  each  phase  of  the  reconstructed 
domain  is  converted  into  an  independent  finite  element  mesh.  The 
diffusional  flow  is  then  computed  for  each  phase  as  a  function  of  its 
bulk  gas  diffusivity  D^k  or  electrical  conductivities  aPulk  (i.e.  D^k 
is  equal  to  the  gas  diffusivity  and  cr|)ulk  is  related  to  electronic  or 
ionic  conductivities).  In  practice,  simulations  are  carried  out  with 
zero  flux  condition  at  the  phase  interface  within  the  domain.  An 
intermediate  boundary  condition  is  applied  onto  the  two  opposite 
faces  along  the  studied  direction  of  the  simulated  volume.  One 
concentration  is  imposed  upon  the  first  face  whereas  the  opposite 
one  undergoes  an  imposed  flux.  The  effective  diffusion  coefficient 
D|gS  (or  the  effective  conductivity  aff)  is  then  deduced  from  the 
numerical  analysis  by  equating  the  macro  homogeneous  flux  to  the 
average  value  computed  from  the  simulation.  The  “apparent”  tor¬ 
tuosity  factor  t  is  then  given  by  the  following  relation  [13]: 

Dgas  =  -^Dgask  anci  aff  =  — °fulk  (i  =  ionic  or  electronic  phase) 

Tgas  Ti 

(1) 

where  e  represents  the  volume  fraction  of  the  analysed  connected 
phase.  It  can  be  mentioned  that  such  definition  of  the  tortuosity 
factor  encompasses  two  distinct  contributions.  The  first  one  is 
related  to  a  pure  “geometrical”  tortuosity.  It  takes  into  account  the 
fact  that  flow  is  hampered  since  the  pathways  for  diffusion  are  not 
straight  but  sinuous  and  intricate.  The  second  contribution  is 
caused  by  some  flux  constriction  effects  in  the  microstructure  since 
cross  section  areas  for  diffusion  are  not  constant  [14-16]. 


where  V(h)  is  the  intersection  volume.  The  differentiation  of  I<(h) 
when  h  =  0  provides  an  exact  determination  of  the  specific  surface 
area,  Sp  [17,18]: 


Sp 


Sarea 

^phase 


with  Sarea 


4  x  map 


h= 0 


(3) 


Such  a  method  offers  the  advantage  to  be  independent  of  the 
reconstructed  phase  digitalization  in  cubic  voxels,  while  a  direct 
computation  of  Sp  is  overestimated  and  depends  of  the  voxelization 

[12]. 

The  covariance  function  C(h)  can  be  used  to  determine  the  mean 
phase  diameter  d.  The  function  represents  the  probability  that  2 
points  of  the  domain  separated  by  a  distance  h  are  included  in  the 
same  phase: 


1  i=N 

C(h)  =P{xe  phase, x+/ig phase}  =  — y^/<(x,y,z)  xk((x+h),y,z) 

i=i 

(4) 

where  N  is  the  number  of  realizations  (i.e.  for  each  distance  h,  N 
corresponds  to  the  number  of  position  x  where  the  relation 
k(x,y,z)  x  k((x  +  h),y,z)  is  calculated).  It  can  be  noticed  that  this  2 
points  probability  function  C(h)  can  be  related  to  the  geometrical 
covariogram  I<(h)  according  to  C(h)  =  K(h)/V(h).  This  function 
tends  to  an  asymptotic  limit  at  a  distance  of  correlation  h*.  It  means 
that  the  2  points  k  and  k(h)  are  not  correlated  anymore.  The  length 
h*  corresponds  to  a  characteristic  dimension  of  the  medium  17,19] 
which  can  ascribed  to  the  mean  particle  diameter  (i.e.  h*  =  d). 
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22.4.  Particle  size  distribution 

Several  methods  can  be  employed  to  determine  the  particle  size 
distribution  20,21].  In  the  present  work,  a  pure  geometrical 
approach  has  been  adopted  [20].  As  illustrated  in  Fig.  2,  the  ana¬ 
lysed  phase  is  represented  by  a  set  of  overlapping  spheres,  each  one 
being  tangent  to  the  complementary  phase.  The  term  p(r2-)  is  then 
defined  as  the  probability  that  a  voxel  belongs  to  the  largest  sphere 
S  of  radius  r2.  In  order  to  compute  this  probability,  voxels  belonging 
to  the  spheres  S(r2)  are  identified  and  their  volumes  are  added 
together  (Fig.  2).  The  resulting  sum  is  normalized  by  the  phase 
volume: 

Pin)  =  T7~  /  X  Vvoxel(X,y,Z)  with/ 

''phase  ^Vphase 

fl  if  (x,y,z)cS(r,-)  f5) 

jo  if  (x,y,z)<£S(rj)  1  ; 

With  Vvoxei  denotes  the  voxel  volume  located  in  (x,y,z)  coordinate. 
The  cumulative  probability,  representing  the  probability  P  that  a 
voxel  belongs  to  a  sphere  of  radius  equal  or  inferior  to  R ,  can  be 
deduced  from  eq.  (5): 


P(0  <  r  <  r„  =  R)  = 


i  =  n  f 

Jfp(ri)  =  /  PSD(r)dr 
1  =  1  i 


(6) 


where  the  Particle  Size  Distribution  (PSD)  corresponds  to  the 
density  probability  function.  This  description  provides  an  exact 
distribution  for  cylinder  or  spherical  morphology.  Nevertheless  it 
has  been  shown  that  it  remains  suitable  for  irregular  medium  [20] 
such  as  the  studied  materials. 


2.2.5.  Triple  phase  boundary  lengths 

The  total  length  of  TPBs  is  calculated  by  counting  all  the  voxel 
edges  in  contact  with  elements  belonging  to  gas,  ionic  and  elec¬ 
tronic  phases.  The  numerical  procedure  implemented  in  Matlab®  is 
the  same  that  the  one  described  in  Vivet  et  al.  [22].  The  computa¬ 
tion  is  performed  according  to  the  three  directions  x,  y  and  z  of  the 
reconstruction.  For  each  direction  (x,  y  or  z-direction),  all  the  voxel 
edges  included  in  the  perpendicular  plane  (yz,  xz  or  xy-plane)  are 
inspected.  If  they  are  identified  as  TPB  segments,  their  lengths  are 
added  to  the  total  length  of  TPB.  Flowever,  such  a  determination 
leads  to  an  over  estimation  since  the  oblique  TPB  line  in  the 
reconstruction  is  constituted  of  orthogonal  segments.  To  take  into 
account  this  shortcoming,  the  corrective  factor  proposed  by  Wilson 
et  al.  [23]  is  applied  by  dividing  the  raw  value  of  TPB  length  by 
1.455.  This  calculated  value  is  referred  as  the  total  TPB  length  since 
it  is  computed  from  the  whole  phases  of  the  reconstructed  volume. 
The  electrochemical  active  length  of  TPB  is  obtained  by  restricting 
the  same  numerical  procedure  to  the  connected  phases. 


2.3.  Electrode  model  description 

The  structural  properties  computed  from  the  3D  reconstructions 
are  introduced  in  an  in-house  electrode  modelling  implemented  in 
Matlab®.  The  model  has  been  already  detailed  in  a  previous  article 
[8].  A  brief  summary  is  given  here  by  highlighting  some  modifi¬ 
cations.  Indeed,  the  model  has  been  adapted  in  order  to  take  into 
account  the  electrode  configuration  made  of  a  Ni-YSZ  composite 
bi-layer  (i.e.  a  functional  layer  layered  on  a  thick  substrate  in  which 
the  electrochemical  reaction  can  potential  extend). 

The  model  is  isothermal  and  considers  one  slice  of  electrode 
embedding  the  functional  layer  and  the  cell  support  (Fig.  3).  Within 


the  functional  layer,  the  electrons  and  ions  fluxes  are  calculated 
according  to  Ohm’s  law  into  the  electronic  and  ionic  conducting 
phases  (cf.  coordinate  system  depicted  in  Fig.  3): 

i'io(z)  =  x  ^  and  ie-  (z)  =  -of  x  (7) 

where  zio  and  ze  are  the  local  electronic  and  ionic  current  densities 
calculated  in  the  Ni  and  YSZ  materials  respectively.  The  terms,  <pio  and 
(pe~  are  related  to  the  electronic  and  ionic  potential  picked  up  in  the 
electronic  and  ionic  conducting  phase.  The  effective  conductivities 
are  expressed  according  to  eq.  (1)  as  a  function  of  the  ‘apparent’ 
tortuosity  factor  and  the  volume  fraction  of  each  phase.  The  bulk 
conductivity  of  the  ionic  8YSZ  phase  is  thermally  activated  and  fol¬ 
lows  an  Arrhenius  law  (aio=8Ysz  (Son-1)  =  466exp(-82, 616/(FT)) 
with  temperature  T  in  I<  [8,24]).  Mass  transfer  of  gas  through  the 
connected  porosity  is  calculated  in  the  frame  of  the  Dusty  Gas  Model 
(DGM)  than  combines  molecular  and  Knudsen  diffusion.  The  effec¬ 
tive  gas  diffusivities  are  still  calculated  by  multiplying  the  bulk  co¬ 
efficients  by  the  ratio  of  the  connected  porosity  to  the  gas  tortuosity 
factor.  Moreover,  it  can  be  mentioned  that  the  Knudsen  coefficient 
depends  on  the  mean  pores  diameter.  In  the  model,  it  is  assumed  that 
the  electrochemical  reactions  can  potentially  occur  in  the  whole 
functional  layer  at  the  active  TPBls.  The  rates  of  the  electrochemical 
reaction  v(z)  are  computed  from  a  Bulter-Volmer  formalism 
considering  a  single  charge  transfer.  The  kinetic  depends  on  the  local 
partial  pressure  of  hydrogen  and  steam,  PH2  (z)  and  PH2o  (z),  computed 
in  the  thickness  of  the  layer: 

v(z)«|koXe  *•  PH2(z)-/<rede  Ph2o(z)|  (8) 

where  R  and  F  denote  the  gas  constant  and  Faraday’s  constant. 
Coefficients  a:ox  and  aTed  correspond  to  symmetry  factors  of  the 
reaction  (aox  +  aired  =  1).  The  function  E(z)  represents  the  elec¬ 
trode  potential  that  is  equal  to  the  difference  between  the  elec¬ 
tronic  and  ionic  potential  taken  in  the  Ni  and  YSZ  phase 
respectively  (i.e.  E(z)  =  <pe\ (z)  -  <pio  (z) ).  Since  the  reaction  takes 
place  at  the  sites  where  ions,  electrons  and  gas  species  meet,  the 
forward  kinetic  constant,  /<ox,  has  been  scaled  with  the  density  of 
‘active’  TPBls,  £TPB1 : 


z=cc- 


Z=FL- 


z=o 


ONi 
O  YSZ 


Fig.  3.  Electrode  geometry  considered  in  the  model.  The  system  of  coordinate  used  for 
the  simulation  is  also  represented. 


398 


F.  Usseglio-Viretta  et  al.  /  Journal  of  Power  Sources  256  (2014)  394-403 


kox  =  Ctpbi  x  /<8x  x  e  rt  (9) 

where  /<gx  denote  the  orientation  factor  and  A Hox  the  activation 
energy  for  charge  transfer.  It  can  be  noticed  that  /<red  is  deduced 
from  the  kox  at  the  Open  Circuit  Potential  (OCP). 

For  the  current  collector  substrate,  it  can  be  observed  in  Fig.  lb 
that  the  structure  of  Ni,  pores  and  YSZ  are  connected  and  contin¬ 
uous  through  the  interface  between  the  functional  layer  and  the 
cell  support.  The  gas,  ions  and  electrons  can  flow  from  one  layer  to 
the  other.  As  a  consequence,  the  electrochemical  reaction  could 
potentially  extend  into  the  bulk  of  the  collector.  In  order  to  estimate 
this  possibility,  the  model  described  in  Ref.  [8]  has  been  modified  in 
such  way  that  the  electrochemical  reaction  could  occur  at  the  active 
TPBls  of  the  substrate.  For  this  purpose,  the  same  set  of  equations 
than  that  used  for  the  functional  layer  are  solved  into  the  cell 
support.  In  this  case,  the  model  parameters  are  expressed  according 
to  the  microstructure  properties  of  the  Ni-YSZ  substrate. 

It  can  be  noticed  that  the  ‘apparent’  electrode  potential  E  cor¬ 
responds  to  the  difference  of  the  electronic  potential  taken  at  the 
free  surface  of  the  electrode  <pel(z  =  CC)  and  the  ionic  potential 
taken  at  the  electrolyte  interface  <?io(z  =  0): 

~E  =  cpel(z  =  CC)  -  cpl0 (z  =  0)  (10) 

The  global  electrode  overpotential  f]  is  then  expressed  as  the 
subtraction  between  E  and  the  OCP  Ei=0:  fj  =  E-EI=0-  This 
expression  of  the  global  electrode  overpotential  can  be  decom¬ 
posed  in  two  distinct  contributions: 

fj  =  ^ohmic  _|_  ^electrochem 

with  ?;ohmic  =  {cpel(z  =  CC)  -  <pel (z  =  z*)}  (11) 

and  jflectrochem  =  {^,(z  =  Z*)  -  cpio(Z  =  0)}  -  Ei=0 

where  z*  is  the  position  in  the  electrode  for  which  the  electro¬ 
chemical  reaction  rates  tends  to  0.  In  other  words,  the  distance  d  = 

| z  |  corresponds  to  the  thickness  of  the  active  layer  since  it  repre¬ 
sents  the  full  extent  of  the  reaction  from  the  electrolyte  interface. 
Therefore,  the  component  ?7ohmic  is  equal  to  a  pure  ohmic  contri¬ 
bution  ascribed  to  the  electrical  resistivity  of  the  inactive  part  of  the 
Ni— YSZ  bilayer;  whereas  the  term  ?|electrochem  represents  the  po¬ 
tential  losses  including  all  the  phenomena  that  hindered  the  re¬ 
action  itself.  It  encompasses  a  combination  of  the  charge  transfer, 
the  gas  transport  up  to  the  TPBls,  and  a  limitation  induced  by  the 
ionic  and  electronic  transport  within  the  active  layer.  This  contri¬ 
bution  corresponds  to  the  so-called  polarisation  resistance  of  the 
electrode  [25]. 

3.  Results  and  data  analysis 

3.1.  Microstructure  quantification 

This  section  of  the  paper  focuses  on  the  determination  of  the 
structural  parameters  from  the  3D  reconstruction  within  cathode 
functional  layer  and  cathode  current  collector.  In  order  to  calculate 
microstructure  properties  with  reasonable  computation  times, 
space  resolution  has  been  downsized  by  increasing  the  voxel  size 
from  25  to  50  nm  (except  for  the  tortuosity  factors  calculated  with  a 
greater  voxel  size  as  further  discussed).  The  procedure  to  decrease 
the  resolution  is  detailed  in  Ref.  [12]. 

The  CFL  reconstructed  volume  lies  just  next  to  the  electrolyte  and 
presents  a  volume  of  27.5  x  27.5  x  10  pm3.  It  has  been  subdivided 
into  9  independent  quasi-cubic  volume  of  9.1  x  9.1  x  10  pm3  taken 
as  Stochastic  Volume  Element  (SVE)  [19,26].  Therefore,  the  CFL 
properties  values  have  been  averaged  on  these  9  realizations. 


Whereas  for  the  CCC,  micro  structural  parameters  have  been  directly 
computed  from  a  single  volume  of  36  x  36  x  36  pm3,  which  is 
slightly  larger  than  the  Representative  Volume  Element  (RVE)  [12]. 


3.1.1.  Volume  fractions,  phase  connectivity  and  TPBls 

Volume  fraction  and  connectivity  for  pore,  YSZ  and  Ni  phase  as 
well  as  TPBls  densities  within  the  CFL  and  the  cathode  substrate  are 
summarized  in  Table  1. 

Difference  in  volume  fractions  between  the  two  layers  is  found 
to  be  significant  as  porosity  increases  from  22.8%  within  CFL  to  46.6 
%  in  the  cathode  substrate.  It  can  be  noticed  that  the  calculated  ratio 
of  solid  volume  fraction  Ni/(Ni  +  YSZ)  (0.63  and  0.59  for  CFL  and 
cathode  substrate,  respectively)  is  in  good  agreement  with  the 
theoretical  value  given  by  the  manufacturer  (0.6  for  50  wt%  Ni  and 
50  wt%  8YSZ).  As  a  consequence,  for  CFL,  9  SVEs  of 
9.1  x  9.1  x  10  pm3  allow  obtaining  good  statistics  (and  ensure  to 
obtain  accurate  mean  values). 

In  Fig.  4,  the  relative  error  in  volume  fractions  calculated  on  the 
CFL  has  been  plotted  as  a  function  of  the  size  of  sub-domains 
extracted  from  the  whole  reconstruction.  It  is  found  that  the  er¬ 
ror  becomes  negligible  for  sub-volumes  larger  than 
14  x  14  x  10  pm3.  This  volume  threshold  defines  the  RVE  size  which 
is  required  to  be  representative  of  the  CFL  fine  microstructure.  As  a 
comparison,  this  volume  is  22  times  lower  than  the  one  requested 
for  the  coarse  microstructure  of  the  cell  support. 

Both  CFL  and  CCC  microstructures  present  a  high  level  of  con¬ 
nectivity.  It  is  greater  than  95%  for  all  phases,  except  for  the  porosity 
within  CFL  which  exhibits  a  lower  value  of  88%.  In  Fig.  5,  the  con¬ 
nectivity  has  been  plotted  for  the  porosity  versus  the  position  in  the 
electrode  thickness.  This  curve  has  been  computed  for  different 
voxel  sizes  ranging  from  50  to  200  nm.  First  of  all,  as  expected,  it 
can  be  noticed  that  the  thickness  of  the  functional  layer  is 
approximately  equal  to  10  pm.  Moreover,  in  the  fine  microstructure 
of  the  functional  layer,  a  voxel  size  of  100  nm,  as  a  maximum,  is 
necessary  to  maintain  the  phase  connectivity.  Indeed,  the  pore 
connectivity  is  found  to  be  strongly  degraded  for  voxel  size  superior 
to  100  nm.  In  the  coarse  microstructure  of  the  CCC,  the  porosity  is 
not  affected  by  the  voxel  size  up  to  200  nm.  These  results  indicate 
that  a  resolution  of  100  nm  is  required  for  the  CFL  while  a 

Table  1 

Volume  fractions  and  connectivity  of  the  phases  plus  total  and  active  TPBls  calcu¬ 
lated  for  CFL,  considering  9  subvolumes  of  9.1  x  9.1  x  10  pm3,  and  cathode  substrate. 
The  standard  deviation  is  calculated  for  CFL  as  well  as  the  extremum  values.  All  TPBls 
values  have  been  corrected  by  the  coefficient  1/1.455. 


CFL 

Cathode 

substrate 

Porosity  volume  fraction  (%)  22.8 

min:  20.3,  max:  26.5,  std.  dev.:  2.1 

46.6 

YSZ  volume  fraction  (%) 

49.0 

min:  46.2,  max:  51.4,  std.  dev.:  1.7 

31.7 

Ni  volume  fraction  (%) 

28.2 

min:  26.7,  max:  29.1,  std.  dev.:  0.8 

21.7 

Porosity  connectivity  (%) 

88.0 

min:  82.7,  max:  93.5,  std.  dev.:  4.3 

99.2 

YSZ  connectivity  (%) 

99.7 

min:  99.6,  max:  99.8,  std.  dev.:  8.2e-2 

99.2 

Ni  connectivity  (%) 

97.1 

min:  95.0,  max:  98.9,  std.  dev.:  1.3 

96.7 

Total  TPBls  (pm  2) 

4.63 

min:  4.42,  max:  4.91,  std.  dev.:  0.16 

3.46 

Active  TPBls  (pm-2) 

3.07 

min:  2.59,  max:  3.42,  std.  dev.:  0.24 

2.61 

Active  TPBls  (pm-2)a 

3.30 

- 

a  Active  TPBls  calculated  on  the  whole  volume  (this  calculation  involves  the 
whole  CFL  volume  and  thus  the  result  is  more  reliable). 
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Fig.  4.  Relative  error  in  volume  fractions  plotted  as  a  function  of  the  domain  size  fi  (for 
an  analysed  volume  of  10  x  £  x  £  pm3  the  constant  dimension  being  along  the  normal 
axis  to  the  interfaces).  Dashed  line  is  the  relative  error  in  the  ratio  of  the  solid  phase 
YSZ/(YSZ  +  Ni).  Volume  analysed  is  extracted  from  the  CFL  reconstruction. 

resolution  of  200  nm  is  sufficient  to  describe  accurately  the  con¬ 
nectivity  of  CCC. 

Total  and  active  densities  of  TPBls  are  both  reported  in  Table  1. 
As  already  mentioned  in  Section  2.2,  the  raw  TPBls  are  corrected  by 
a  coefficient  of  1/1.455  because  of  the  discretisation  in  cubic  voxels. 
As  expected,  the  CFL  shows  higher  active  TPBls  density  than  CCC 
(i.e.  3.3  pm-2  for  CFL  and  2.6  pm-2  for  CCC).  It  can  be  remarked  that 
the  value  of  total  TPBls  within  CFL  (i.e.  4.63  pm-2)  is  very  close  to 
the  one  given  by  Joos  et  al.  [27  on  the  same  kind  of  electrode  (6.79/ 
1.455=4.66  pm-2). 

3.1.2.  Specific  surface  area,  mean  particle  size  and  particle  size 
distribution 

Specific  surface  area,  mean  particle  diameters  for  pore,  YSZ  and 
Ni  phases  within  CFL  and  cathode  substrate  are  summarized  in 


Position  along  the  normal  axis  with  the  interfaces  (pm) 

Fig.  5.  Porosity  connectivity  calculated  on  the  whole  volume  and  here  plotted  on  each 
slice  along  the  normal  axis  with  the  interfaces.  Data  are  taken  from  the  CFL  recon¬ 
struction  with  a  part  of  the  electrolyte  on  the  left.  Data  has  been  plotted  with  different 
voxel  size  ranging  from  50  to  200  nm. 


Table  2.  Note  that  these  values  have  been  computed  with  geometric 
covariograms  and  covariance  functions  (see  Section  2.2). 

CFL  systematically  exhibits  smaller  particles  size  compared  to 
CCC,  as  deduced  form  the  covariance  functions  plotted  in  Fig.  6.  This 
is  especially  true  for  the  porosity  with  a  mean  particle  size  of 
0.94  pm  in  the  functional  layer  versus  2.6  pm  in  the  support.  Note 
that  the  large  difference  in  porosity  size  between  CFL  and  CCC  can 
be  clearly  observed  on  the  3D  reconstruction  (Fig.  1).  The  presence 
of  large  pores  in  the  substrate  after  sintering  results  in  the  use  of 
pore-former  agents  during  the  manufacturing  process  [28,29]. 

Specific  surface  area  for  porosity  is  more  than  twice  within  CFL 
compared  to  CCC  (i.e.  9.09  versus  4.24  pm-1).  This  result  is  in  good 
agreement  with  the  difference  in  pores  diameters  between  the  two 
layers. 

Aside  from  the  method  based  on  geometric  covariogram  (see 
Section  2.2),  specific  surface  area  has  been  also  computed  from  a 
direct  analysis  of  3D  reconstruction.  In  this  case,  the  specific  surface 
area  is  calculated  by  adding  all  the  voxel  faces,  which  are  located  at 
the  boundary  between  the  phase  and  the  complementary  volume 
[12].  However,  this  method  requires  a  corrective  factor  to  obtain  the 
real  surface  since  the  structure  is  described  by  cubic  elements.  A 
factor  of  2/3  has  been  considered  in  this  work  that  is  rigorously 
applicable  only  for  a  sphere  described  by  cubic  voxels.  In  spite  of 
this  simplification,  specific  surface  areas  computed  from  the  two 
methods  are  found  to  be  in  good  agreement.  Indeed,  relative  error 
between  the  two  methods  is  less  than  7.5%.  As  a  consequence,  the 
direct  method  based  on  the  voxel  faces  summation  provides  good 
approximation  for  SOEC  electrode  microstructure. 

PSDs  for  CFL  and  CCC  are  provided  in  Fig.  7.  It  has  to  be  remarked 
that  PSD  algorithm  gives  smaller  mean  particle  diameters  than  the 
covariance  functions.  Indeed,  the  continuous  PSD  algorithm  used  in 
this  work  leads  to  underestimate  the  real  mean  particle  diameter  as 
already  discussed  by  Holzer  et  al.  [30]. 

Whatever  the  electrode  layer,  histograms  of  solid  phases  are 
quite  similar.  They  present  a  unimodal  distribution  with  a  relative 
narrow  peak.  Conversely,  two  distinct  histograms  are  found  for 
porosity  of  CFL  and  CCC.  The  cathode  support  exhibits  a  wide  dis¬ 
tribution  of  diameters  while  the  functional  layer  is  centred  on  a 
narrow  peak.  As  discussed  previously,  the  wide  distribution  of 
pores  size  in  the  substrate  could  be  explained  by  the  use  of  pore- 
former  agents,  which  should  result  in  large  voids  after  sintering. 

3.1.3.  Tortuosity  factor 

A  sensitivity  analysis  has  been  performed  on  the  tortuosity 
factor  by  changing  the  voxel  size.  For  the  functional  layer,  it  has 
been  found  that  the  size  of  the  discretization  must  be  lower  than 
<100  nm  to  compute  accurate  tortuosities.  For  the  cathode  support, 
the  same  condition  on  voxel  size  is  requested  for  the  Ni  and  YSZ 
phases  (<100  nm).  Regarding  the  gas  phase  of  the  substrate,  it  has 
been  found  that  voxels  of  200  nm  are  small  enough  to  compute 
stable  tortuosity  factors.  It  can  be  remarked  that  these  statements 
are  consistent  with  the  sensitivity  analysis  performed  on  phase 
connectivity.  Moreover,  a  previous  study  [12]  performed  on  a 
typical  Ni— YSZ  support  has  already  shown  that  a  discretisation  of 
240  nm  is  sufficient  to  determine  the  gas  tortuosity  factor  (with  an 
error  lower  than  2.5%).  Explanation  for  the  dependence  of  tortu¬ 
osity  factor  on  voxel  size  derives  from  the  particle  size  distributions 
between  the  phases.  Porosity  of  the  cathode  substrate  exhibits  a 
non-homogeneous  particle  size  distribution  (Fig.  7),  that  can  be 
simplified  with  a  fine  and  a  coarse  pore  network.  Gas  species 
diffuse  mainly  through  the  coarse  network  of  the  structure  while 
the  fine  network  has  only  a  limited  contribution.  The  increase  in 
voxel  size  leads  to  the  removing  of  a  part  of  the  fine  network 
(Fig.  8),  resulting  to  a  slight  increase  in  the  tortuosity  factor. 
Conversely,  the  functional  layer  and  the  solid  phases  of  the  support 
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Table  2 

Specific  surface  area  Sp  and  mean  particle  diameter  dmean  of  the  phases  calculated  for 
CFL,  considering  9  subvolumes  of  9.1  x  9.1  x  10  pm3,  and  cathode  substrate.  The 
standard  deviation  is  calculated  for  CFL  as  well  as  the  extremum  values.  The  specific 
surface  area  has  been  calculated  with  the  geometric  covariogram  and  the  mean 
particle  diameter  with  the  covariance  function. 


CFL 

Cathode 

substrate 

Sp  Porosity  (Fm  ) 

9.09 

min:  7.88,  max:  10.19,  std.  dev.:  0.81 

4.24 

Sp  ysz  (nm-1) 

5.50 

min:  5.32,  max:  5.69,  std.  dev.:  0.12 

6.01 

Sp  Ni  (|im-1) 

6.54 

min:  6.30,  max:  6.69,  std.  dev.:  0.14 

5.87 

^mean  Porosity  (|-im) 

0.94 

min:  0.65,  max:  1.7,  std.  dev.:  0.35 

2.6 

dmean  YSZ  (hm) 

0.73 

min:  0.7,  max:  0.75,  std.  dev.:  0.03 

1.1 

dmean  Ni  (hm) 

0.83 

min:  0.75,  max:  0.9,  std.  dev.:  0.06 

1.2 

except  for  the  porosity  of  the  CCC  for  which  the  resolution  has  been 
taken  to  200  nm.  The  obtained  values  along  the  perpendicular  di¬ 
rection  to  the  electrolyte  (i.e.  z-direction:  cf.  Fig.  1)  are  reported  in 
Table  3.  Ni  and  YSZ  structure  present  lower  tortuosity  factor  within 
CFL  than  substrate.  Since  connectivity  and  PSD  are  roughly  identical 
between  the  two  layers,  that  effect  could  be  ascribed  to  the  dif¬ 
ference  in  volume  fraction  between  the  two  layers. 

Unlike  the  solid  phases,  porosity  exhibits  a  tremendous  increase 
in  tortuosity  factor  from  2.22  for  the  substrate  to  10.15  for  the 
functional  layer.  These  differences  can  be  explained  by  both  the 
change  in  volume  fractions  and  pore  size  distribution. 

3.2.  Ni-YSZ  electrode  modelling 

Electrode  simulations  have  been  carried  out  with  the  micro- 
structural  properties  determined  in  the  previous  section  for  the  Ni- 
8YSZ  bi-layer.  The  kinetic  parameters  introduced  in  the  model  are 
provided  in  Table  4. 


present  only  a  fine  microstructure  (Fig.  7)  that  has  to  be  described 
by  a  higher  space  resolution. 

Since  homogenization  computations  are  time  consuming,  the 
tortuosity  factors  have  been  calculated  with  a  resolution  of  100  nm, 


Fig.  6.  Covariance  function  plotted  as  a  function  of  the  distance  h  calculated  (a)  on  CFL 
and  (b)  on  the  cathode  substrate.  For  the  sake  of  visibility  the  covariance  is  plotted  only 
for  one  subvolume  of  9.1  x  9.1  x  10  pm3  for  CFL. 


3.2.  t.  Orientation  factor  kgx 

For  the  best  knowledge  of  the  authors,  there  is  no  value  reported 
in  the  open  literature  regarding  to  the  kinetic  parameter  kgx  related 
to  a  Ni-8YSZ  cathode  operating  in  steam  electrolysis.  As  a  conse¬ 
quence,  the  pre-exponential  factor  kgx  has  been  determined  in  the 
present  study  so  that  the  ‘apparent’  exchange  current  density  of  the 
electrode  is  equal  to  i0  =  537  mA  cm-2  at  T  =  800  °C.  Indeed,  this 
value  has  been  experimentally  determined  in  another  study  for  the 
same  cell  that  the  one  studied  in  this  paper  [32].  It  is  reminded  that 
the  parameter  io  is  associated  to  the  electrode  activation  over¬ 
potential,  which  is  expressed  as  follows  [5,32]: 

^act  =  Tsinh~m  (12) 

The  overpotential  7]act  includes  all  the  phenomena  that  hinder  the 
electrochemical  process  except  the  limitation  due  to  the  gas 
transport  (that  is  ascribed  to  the  concentration  overpotential 
[31,32]).  The  electrode  overpotential  ^eiectrochem^  which  can  be 
calculated  with  the  electrode  model  (see  eq.  (11)),  is  then  equal  to 
the  activation  overpotential  ?7act  only  at  low  current  density  for 
which  the  gas  diffusion  is  not  a  limiting  step  [5,33].  As  a  conse¬ 
quence,  the  methodology  used  to  estimate  /<gx  has  been  the 
following:  a  set  of  simulations  have  been  carried  out  at  T=  800  °C 
and  at  low  current  density  (i.e.  | ie  |  <  0.5  A  cm-2).  The  simulations 
were  carried  out  by  changing  kgx  up  to  obtain  a  satisfactory 
adjustment  between  the  electrode  overpotentials  ^electrochem  and 
the  activation  ones  (i.e.  ^act  calculated  with  eq.  (12)  by  considering 
an  ‘apparent’  exchange  current  density  of  537  mA  cm-2  32]).  As 
shown  in  Fig.  9,  the  best  fit  is  obtained  for  /<qx  =  29  mol  s-1  m-1.  It 
can  be  mentioned  that,  even  at  the  highest  investigated  current 
density,  the  pressure  drop  across  the  Ni-8YSZ  bilayer  remains  fairly 
low  meaning  that  the  concentration  overpotentials  can  be  neglec¬ 
ted  (i.e.  APh2o  =  0.02  at  ie-  =  -0.5  A  cm2).  This  statement  is 
consistent  with  previous  studies  [5,33]  that  have  shown  that  sig¬ 
nificant  concentration  overpotentials  arise  for  typical  electrode 
supported  cells  only  at  high  current  density  (|ie-|  >  1  A  cm-2). 

3.2.2.  Electrode  activation  energy 

Since  the  electrochemical  processes  are  thermally  activated,  the 
‘apparent’  exchange  current  density  is  usually  expressed  according 
to  an  Arrhenius  law  [5]: 

io(T)  =  i00  x  e-&  (13) 

where  Ea  is  the  electrode  activation  energy  which  is  classically 
employed  in  a  ‘macro’  model  developed  at  the  cell  or  stack  level 
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Fig.  7.  Particle  size  distribution  (a)  on  CFL  and  (b)  on  the  cathode  substrate.  Voxel  size  is  reduced  to  50  x  50  x  50  nm3  in  order  to  describe  the  smallest  spheres  within  CFL  whereas 
voxel  size  is  increasing  to  100  and  200  nm  respectively  for  solid  phases  and  porosity  within  cathode  substrate  to  get  reasonable  time  calculations.  The  dashed  lines  indicate  the 
mean  values  of  the  given  distribution. 


[5,32-34].  Nowadays,  despite  its  great  importance,  no  specific 
study  has  been  devoted  to  determine  this  parameter  for  a  Ni— YSZ 
cathode. 

In  order  to  estimate  a  reliable  value  for  Ea,  the  electrode  over¬ 
potential  ^eiectrochem  jias  been  computed  with  the  cathode  model  at 
T  =  750  °C  and  T  =  850  °C.  As  previously  mentioned,  the  calcula¬ 
tions  are  restricted  to  low  current  cases  to  avoid  any  limitation  due 
to  the  gas  transfer.  The  simulated  cathode  overpotential  has  been 
then  fitted  by  eqs.  (11)  and  (12)  by  considering  an  exchange  current 
density  io  of  537  mA  cm-2  at  T  =  800  °C.  The  fit  result  is  given  in 
Fig.  9.  It  is  found  that  the  dependence  on  temperature  of  the 
electrode  overpotentials  is  accurately  described  by  relations  (12) 
and  (13)  for  an  activation  energy  of  Ea  =  120.5  kj  mol-1.  This 
value  encompasses  the  activations  energies  associated  to  the 
thermally  activated  phenomena  controlling  the  whole  electro¬ 
chemical  process.  It  can  be  noticed  that  the  value  obtained  in  this 
study  is  in  good  agreement  with  the  reported  activation  energies  by 
Costamagna  and  Honegger  35]  (i.e.  Ea  =  120  and  140  kj  mol-1  for 
an  Ni-8YSZ  composite  operated  as  SOFC  anode  under  hydrogen  at 
800  °C). 

At  low  current  density,  it  can  be  reasonably  expected  that  the 
electrode  response  is  mainly  governed  by  the  kinetic  of  the  elec¬ 
trochemical  reaction  and  the  oxygen  ions  transport  in  the  YSZ 
structure  of  the  active  layer  [8].  Indeed,  the  contribution  of  the 
species  depletion  along  the  electrode,  or  the  concentration  over¬ 
potential,  can  be  neglected  in  this  case.  However,  this  work  shows 
that  the  electrode  activation  energy  is  not  a  single  combination  of 
the  activation  energies  related  to  the  ionic  conduction  (i.e. 
83  kj  mol-1)  and  charge  transfer  on  Ni-8YSZ  material  (i.e. 
201  kj  mol-1).  Thanks  to  a  micro  model,  only  the  correct  descrip¬ 
tion  of  the  reaction  delocalization  through  the  functional  layer  and 
its  microstructure  properties  is  able  to  supply  reliable  values  of  the 
effective  exchange  current  density. 


3.2.3.  Thickness  of  the  active  layer  d  =  |z*  | 

Thanks  to  the  present  study,  the  microstructure  properties  of 
electrode  and  its  kinetic  parameters  are  known.  As  a  consequence, 
it  becomes  possible  to  assess  with  the  electrode  model  the  real 
extent  of  the  electrochemical  reaction  into  the  Ni-8YSZ  bilayer.  For 
this  purpose,  the  local  volume  specific  current  j(z)  has  been 
computed  versus  the  position  z  in  the  electrode  (Fig.  10).  It  is 
reminded  that  j(z)  represents  the  electrochemical  source  terms  that 
are  produced  in  the  active  layer.  Indeed,  the  local  specific  current 
j(z)  is  linked  to  the  kinetic  rates  of  the  electrochemical  reaction  v(z) 
according  to  j(z)  =  2 F  x  v(z).  As  shown  in  Fig.  10,  it  is  found  that 
source  terms  j(z)  associated  to  the  electrochemical  reactions  occur 
in  the  whole  functional  layer;  whereas  the  extent  of  the  reactions 
into  the  substrate  is  very  limited.  This  result  can  be  explained  by 
two  main  reasons: 

(i)  Firstly,  the  decrease  of  the  electrochemical  reaction  rates  is 
very  steep  in  the  depth  of  the  functional  layer  (Fig.  10).  At 
10  pm  from  the  electrolyte  interface,  the  electrochemical 
reaction  rates  start  to  be  limited  (they  represent  only  ~  5%  of 
the  maximum  value  obtained  at  the  interface  with  the 
electrolyte  (z  =  0)).  This  sharp  decrease  of  the  electro¬ 
chemical  reaction  is  mainly  explained  by  the  difficulty  of 
oxygen  ions  to  diffuse  over  a  large  distance  into  the  func¬ 
tional  layer  [8].  It  means  that,  for  the  studied  microstructure, 
the  delocalisation  of  the  electrochemical  reaction  beyond  an 
active  layer  of  10  pm  is  very  limited. 

(ii)  The  previous  phenomenon  is  reinforced  by  the  poor  elec¬ 
trochemical  activity  of  the  cell  substrate  because  of  its  coarse 
microstructure.  Indeed,  the  concentration  of  active  sites  is 
significantly  lower  in  the  cell  support  in  comparison  to  the 
functional  layer  (i.e.  ^^pP0,t/X™cti°na| 'a^r  =  0  8)  More_ 
over,  the  oxygen  ions  mobility  through  the  YSZ  structure  of 


402 


F.  Usseglio-Viretta  et  al.  /  Journal  of  Power  Sources  256  (2014)  394-403 


Fig.  8.  Bidimensional  view  extracted  from  the  middle  of  the  cathode  substrate  reconstruction.  White  is  the  percolating  porosity  phase  and  black  is  the  complementary  volume. 
Voxel  size  is  (a)  100  x  100  x  100  nm3,  (b)  200  x  200  x  200  nm3  and  (c)  400  x  400  x  400  nm3.  The  fine  percolating  network  is  preserved  with  a  discretization  up  to  a  voxel  size  of 
200  nm  whereas  it  has  almost  been  removed  for  a  voxel  size  of  400  nm. 


Table  3 

Tortuosity  factor  calculated  for  CFL,  considering  9  subvolumes  of  9.1  x  9.1  x  10  pm3, 
and  cathode  substrate.  The  standard  deviation  is  calculated  for  CFL  as  well  as  the 
extremum  values.  The  tortuosity  factor  has  been  calculated  along  the  axis  normal  to 
the  interfaces.  Voxel  size  is  200  nm  for  porosity  within  cathode  substrate  and 
100  nm  for  all  others  phases. 


CFL 

Cathode  substrate 

Porosity 

10.15 

min:  7.26,  max:  15.06,  std.  dev.:  2.88 

2.22 

YSZ 

2.28 

min:  1.97,  max:  2.88,  std.  dev.:  0.33 

3.63 

Ni 

5.50 

min:  4.40,  max:  7.68,  std.  dev.:  0.99 

8.65 

Table  4 

Kinetic  parameters  for  Ni— YSZ  electrode  simulations. 


Electrode 

Activation  energy 

Orientation 

Symmetry 

material 

for  charge  transfer  A Hox 
a  mol-') 

factor  k°Qx 
(mol  s~H  m"1) 

factor  ac  and  aa 

Ni-8YSZ 

201  x  103  [31] 

29a 

0.5 

a  Value  determined  in  this  study. 


Fig.  9.  Electrode  overpotential  7yelectrochem  simulated  for  the  Ni-8YSZ  cathode  bi-layer 
in  10%  H2  and  90%  H20.  The  microstructure  properties  taken  for  the  simulations  are 
provided  in  Section  3.1.  The  simulated  points  are  compared  to  the  activation  over¬ 
potential  ?7act. 
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Fig.  10.  Delocalization  of  the  electrochemical  reaction  for  the  Ni-8YSZ  cathode  bi-layer 
in  10%  H2  and  90%  H20  (at  T  =  800  °C  for  an  electrode  overpotential  of  0.02  V).  The 
microstructure  properties  taken  for  the  simulations  are  provided  in  Section  3.1. 


the  cell  support  is  very  limited.  Indeed,  the  tortuosity  factor 
of  this  ionic  conducting  phase  in  the  current  collector  sub¬ 
strate  is  much  higher  than  the  one  related  to  the  functional 
layer  (i.e.  TYS2Support/TYSzCti°nal  layer  -  1.6). 

Therefore,  it  can  be  inferred  from  this  analysis  that  the  delo¬ 
calization  of  the  electrochemical  reaction  within  the  cell  support  is 
very  low  meaning  that  its  participation  to  the  electrochemical  re¬ 
actions  can  be  neglected.  As  a  consequence,  it  can  be  assumed  that 
the  dimension  of  the  active  layer  can  be  approximately  restricted  to 
the  functional  layer  (i.e.  d  =  \z  =  FL|). 

4.  Conclusion 

Microstructure  of  a  cermet  Ni-YSZ  bi-layer  has  been  analysed 
on  the  basis  of  3D  reconstructions  obtained  on  both  functional  layer 
and  cell  support.  Thanks  to  the  large  field  of  view  of  the  re¬ 
constructions  ( —  50  pm),  microstructure  properties  have  been 
quantified  with  a  high  level  of  confidence.  Indeed,  it  has  been 
shown  that  the  analysed  volumes  are  statistically  representative  of 
the  electrode  bi-layer  microstructure. 

For  the  first  time,  all  the  microstructure  properties  of  gas,  ionic 
and  electronic  phases  have  been  determined  for  the  studied  typical 
Ni-YSZ  cermet.  Connectivity  of  each  phase,  mean  particles  di¬ 
ameters,  specific  surface  area,  particle  size  distributions,  tortuosity 
factors  and  densities  of  TPBls  have  been  quantified  for  both  func¬ 
tional  layer  and  cell  support. 

Microstructure  parameters  have  been  introduced  in  an  SOEC 
cathode  micro  model  that  takes  into  account  the  specific  configu¬ 
ration  of  the  Ni-YSZ  composite  double-layer.  Simulations  have 
shown  that  the  extent  of  the  electrochemical  reaction  in  the  sup¬ 
port  is  very  limited  and  can  be  neglected.  The  SOEC  cathode  ex¬ 
change  current  density  has  been  computed  as  a  function  of 
temperature.  It  has  been  found  that  the  combination  of  ionic  con¬ 
duction  and  charge  transfer  all  along  the  volume  of  the  electrode 
results  in  an  apparent  activation  energy  of  120.5  kj  mol-1. 
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